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ABSTRACT 


Experimental work was performed to study the convective heat transfer 
characteristics from a cyhnder in a strong zero-mean oscillatory flow represented by 
an acoustic field. Two different flow regimes are discussed; that in which laminar, 
attached flow around the cylinder is present, and that in which instabilities, such as 
vortex shedding occur. The experiment utilizes a steady state measurement mediod. 
A transition from the laminar to the unstable regime was observed to occur at a 
streaming Reynolds number of approximately 240. Within the laminar regime, tiie 
transition from “intermediate” to “large” values of the streaming Reynolds number 
occurs at approximately 130. Heat transfer results for large values of the streaming 
Reynolds number in the laminar regime closely match the present theory (less than 
13% error). Correlations were developed to relate the heat transfer rate to die 
streaming Reynolds number in the unstable regime. This work would fmd apphcation 
in the design of heat exchangers for a thermoacoustic engine. 
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I. INTRODUCTION 


The science of convective heat transfer in an acoustic field, while still in its infancy, 
presents many new and exciting possibilities for application in the future. Different 
experiments and theory have proven that under the correct circumstances there are certain 
desirable heat transfer characteristics affecting an object which is immersed in a strong 
acoustic field. A complete analysis concerning the processes and effects of the related heat 
transfer phenomena, though, is lacking and desperately needed. 

Although thermoacoustics have already been applied to some advanced heat transfer 
designs, for instance, the thermoacoustic space refrigerator and other thermoacoustic 
cryocoolers developed at the Naval Postgraduate School, there has yet to be developed 
anything that can compete on an economic level with what is currently marketed today. The 
efficiencies obtained so far have been quite low, requiring nearly twice the power of a 
conventional vapor compression refrigerator. When the fundamentals behind the 
thermoacoustic phenomenon and the related heat transfer characteristics are completely 
understood, breakthroughs can occur which could allow industry to move ahead and apply 
these techniques on an every day basis toward a variety of common uses. 

To properly model and control the parameters which impact upon the heat transfer 
behavior in a thermoacoustic engine, it would be advantageous if the various flow regimes 
(e.g., turbulent vs. laminar) in the engine could be isolated and analyzed in detail. Further 
information in this regard may be obtained through a parametric analysis of a suitable model 
problem by which a measure of the importance (or rather a magnitude of the effect each 
parameter has on the process as a whole) can be determined. This is an important element 
of the modeling process and requires study since by understanding the impact that each 
individual parameter makes upon the thermoacoustic process as a whole, it may be possible 
to predict the changes in the heat transfer characteristics as individual components are varied. 

The work contained within provides an experimental study of some of the dominant 
heat transfer properties of a particular model problem that may be encountered in 
thermoacoustic engines. The model problem chosen is one of convective heat transfer from 
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a cylinder in a zero-mean oscillatory flow. The flow is representative of the acoustic 
standing wave in a thermoacoustic engine whereas the cylinder represents a tube or other 
component that may be present in such an engine. 

The work involves a correlation of experimental heat transfer data in terms of a 
suitable Nusselt number (Nu) with other appropriate dimensionless parameters in the 
problem, such as the streaming Reynolds number (R^), which itself is a function of length 
scales, pressure ratios and frequency parameters and, of course, of the Prandtl number (Pr). 
Through use of high power standing resonant acoustic waves in a cylindrical chamber, a 
high-intensity internal oscillatory flow is established. Under these conditions, the heat 
removal rate from a thin cylindrical heating element immersed in the acoustic signal will 
supply data necessary to arrive at some basic conclusions as to heat transfer phenomena 
occurring around a cylinder. This work focuses upon two different flow regimes; one in 
which laminar, attached flow around the cylinder at large values of the streaming Reynolds 
number is present, and the second in which vortex shedding and other instabilities in the flow 
are expected to occur at the cylinder surface. The resultant experimental data additionally 
provides guidelines for determining when the flow transitions from one regime to the other. 
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II. BACKGROUND 


A. HISTORICAL 

It has long been understood that a large temperature gradient along the length of a 
cylindrical tube can, under certain suitable circumstances, spontaneously excite the fluid into 
oscillations strong enough to create audible sound. Glass blowers provided the earliest 
accounts of this acoustic effect. They found that when one end of a glass tube was placed 
in a furnace, the temperature difference between the end of the tube in the furnace and the 
end still under ambient conditions created an audible tone which was emitted from the open 
end of the tube. Even though early scientists knew of this effect, it was merely considered 
to be more of an oddity than a scientific discovery which might have useful implications. 
In fact, the earliest accounting of what may have been the first ever thermoacoustic engine, 
the Sondhauss tube (Figure 1), was described by Sondhauss (1850) himself as the “glowing 
glass harmonica”. Lord Rayleigh (1945) gave the first good qualitative analysis of the 
Sondhauss tube in 1896 in which he described the mechanism behind the effect and posed 
the theory that mechanical work could be obtained from the vibrations, or oscillations, which 
were being created by the temperature gradient along the tube. 



3 





An extension of research into Sondhauss tubes was conducted by Carter (1988) and 
included placing a stack of plates at an appropriate point in the tube. The plates included hot 
heat exchanger strips at one end and cold heat exchanger strips at the other. These improved 
the effect of the Sondhauss tube, and inspired another scientist, Feldman (1966) to conduct 
similar research which consequently resulted in an oscillator which produced 27 W of 
acoustic power from 600 W of thermal power. All of this occurred in the 1960s. 

It wasn’t until much later that it was postulated that the effect could somehow be 
reversed, i.e., a temperature gradient could be created along a tube using a poweiM, 
resonating acoustic signal as the driving force. It is only more recently that advances in 
acoustic technology have allowed serious research toward this goal. 

B. RECENT WORK 

The heat transfer effect related to this thermoacoustic phenomenon has since be^ 
used in the creation of heat pumps and has been explored by several scientists. The results 
of these attempts include what was termed a pulse-tube refrigerator as developed by Gifford 
and Longsworth (1966). This refrigerator utilized low frequency, high amplitude oscillations 
to excite the gas in a tube and create a cooling effect along the surface. The invention of 
“modem” thermoacoustic refrigeration occurred in the early 1980s at Los Alamos National 
Laboratory. It was in essence a modification to the work that Carter did, using stacks of 
plates with a much smaller temperature gradient. Additional engineering developments by 
others, such as the work at the Naval Postgraduate School with a thermoacoustic reftigeratOT 
intended for use on the space shuttle, led to increases in efficiency, as well as to an increase 
in commercial interest and development. Currently, there are major projects ongoing in 
several countries as outlined by Swift (1995), including a prototype food refrigerator based 
upon the Naval Postgraduate School’s work being built in the Republic of South Africa. The 
Ford Motor Company has developed its own version of a thermoacoustic refrigerator while 
the Tektronix Corporation is working towards a pulse-tube type of refrigerator to be used for 
cooling electronics to cryogenic temperatures. 

All of these attempts are especially significant given the current stigma surrounding 
the environmentally hazardous use of CFC’s. Even though thermoacoustic refrigeration is 
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still not as efficient as that of current energy efficient vapor compression models, there is a 
growing demand for something to take their place. By advancing our knowledge-base in this 
area, and further incorporating a new understanding of how increases to the efficiency of 
thermoacoustic designs can be made, that it may well be possible that a new thermoacoustic 
revolution is in our future. 

C. THERMOACOUSTIC PROCESS 

The basic process behind a thermoacoustic engine can be best described by the model 
in Figure 2. The upper portion of the figure shows a sound chamber with an acoustic driver 
at the left end which is used to create a resonant, standing wave in the chamber. At an 
appropriate point within the chamber, a thermoacoustic stack (Figure 3) is placed with a heat 
exchanger (Figure 4) on either side of it. The flow of heat is from right to left in the figure. 
The process of heat transport across the plate is illustrated in the lower portion of the figure. 
The fluid within the chamber will oscillate due to the acoustic wave, traveling from a point 
of low pressure to one of high pressure, gaining and losing energy during each half cycle. 

For instance, a parcel of gas at temperature Tq at low pressure moves along line 1, 
increasing in pressure and temperature until it reaches T-h- (has gained two units of heat). 
At that point, it loses one unit of heat to the plate, thereby reducing its temperature to T+. 
As the parcel of gas continues through the second half of the oscillatory cycle, it decreases 
in pressure, and it loses two unit of heat, dropping to T-. It is then able to retrieve a single 
unit of heat from either the plate at the right hand side of the cycle, or from a heat exchanger 
at that end. In effect, a “bucket-brigade” of little parcels of gas is formed as the heat 
transport mechanism. It is the heat exchangers on either end of the thermoacoustic stack and 
the heat transport mechanisms involved with them that this experiment intends to analyze. 

D. HEAT TRANSFER IN ACOUSTIC FIELDS 

The early 1950's saw an increased interest towards understanding the heat transfer 
behavior in oscillatory flows and an earnest effort towards understanding the possible 
benefits thereof began. Richardson (1967) produced the first significant contribution to this 
field by providing a coherent and detailed account on the general nature of heat transfer in 
oscillatory flows. He gave concise documentation on how sound and vibration fields had 
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made an impact in the heat transfer field up to that point. He additionally conducted a 
preliminary study on how the effect can be related to that of a more traditional fluid- 
mechanical analysis. Davidson (1974) further expanded upon Richardson’s work, analyzing 
the heat transfer behavior of cylinders in oscillatory flow. That occurred in the early 1970's, 
and since then, very little work on the subject can be found in literature. 

Most recently, though, Mozurkewich (1995) performed heat transfer experiments in 
an acoustic field utilizing a transient analysis of a heated wire. His results, although 
informative in some respects, were lacking in data within the most basic flow regimes and 
his conclusions left some doubt to the reader. 

When first attempting to analyze the process behind convective heat transfer in an 
acoustic field, it is often simpler to think of the heat transfer phenomenon as something akin 
to that of forced convection due to a steady mean flow. Initially, that there is a separate 
power source placed away from the test object which produces a disturbance in die fluid 
medium in which that object is immersed. In forced convection, that power source may be 
considered to be a fan or a pump which creates a pressure gradient, which in turn causes a 
steady flow of fluid. In the current problem, the power source is instead an acoustic driver 
which causes an oscillatory (or vibrational) type of time-periodic flow around the 
being considered. This oscillatory flow has a zero-mean and results in no net through flow. 
Before analyzing this flow for heat transfer characteristics, though, it is necessary to first iK)te 
which aspects of the acoustic field dominate. 

A resonant, standing wave acoustic field excited across the ends of a closed 
cylindrical chamber has very distinct properties. Of particular interest is the fact that when 
such a field exists, the pressure along the length of the chamber varies sinusoidally sudi that 
a point of maximum pressure occurs at the rigid end termination at the opposite end of die 
chamber from which the acoustic signal is being generated. In addition, depending uptm die 
frequency being used, there may be more than one zero-crossing, or pressure minimum, 
along the length of the chamber (Figure 5). 
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At resonance, the acoustic velocity is out of phase with the pressure. For instance, 
at a point of minimum pressure, a pressure node, the particle oscillations will be at their point 
of maximum velocity, a velocity antinode. The reverse then also holds tme, (i.e., a pressure 


Circular Sound 



Velocity Profile 
Pressure Profile 


Figure 5. Pressure and Velocity Profiles in a Circular Tube in the 
Presence of an Acoustic Signal. 


antinode can also be designated a velocity node). 

To understand the reason why this relationship is extremely valuable to this research, 
it is useful to again refer to the forced convection model. In general, as the fluid velocity in 
a forced convection application increases within a particular regime, so does the heat transfer 
coefficient. The case is essentially the same for heat transfer in an acoustic field, and hence 
it is expected that the heat transfer rate would also increase as the particle velocity increases. 
Thus, for experimental purposes it is of prime importance that the test object be placed at a 
velocity anti-node to get the maximum effectiveness out of the process for a given acoustic 
signal. However, in contrast to the forced convective mean flow case, the current issue with 
oscillatory flow is considerably more complicated due to the wide range of flow parameters, 
and hence flow patterns and heat transport mechanisms, that can result. 
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II. EXPERIMENT 


A. INTENT 

As with any other advance in technology or new scientific discovery with which 
engineers desire to predict and quantify results in some manner, it is best to begin by first 
breaking the process down into its basic component parts. By completing an analysis for the 
simplest version of the model in question, a stepping stone will be established upon which 
the analysis of more complicated scenarios can be built. It stands to reason that this is how 
the analysis of the thermoacoustic heat transfer process should also begin. 

For an initial starting point, the analogy once again to forced convection heat transfer 
is used. The problem of an isolated cylinder in a mean cross flow is well documented and 
understood, and has in turn been used to develop correlations for flow over a collection of 
cylinders, such as a tube bank, a more practical application as is evident from any basic heat 
transfer textbook. This is the motivation for a study of the behavior of simple shapes such 
as a cylinder placed in an acoustic field. It is hoped that with this knowledge for an isolated 
cylinder, it would be possible to extend the solution to other models. 

It is only necessary then to concentrate on the evaluation of the acoustic signal itself 
in terms of its many different parameters. But before this analysis may begin, the type of 
acoustic flow around the cylinder must be established. To allow corroboration with 
established theory, the flow pattern initially desired is that of basic laminar, incompressible 
flow where well understood streaming patterns are the principal forms of flow and heat 
transport. This meets the requirement for maintaining the most simplified version of flow 
for the analysis. 

The geometry of how the test cylinder is placed with reference to the acoustic signal 
is also of utmost importance. A theoretically perfect scenario would have the test cylinder 
situated normal to a unidirectional sound field with no interference from the surroundings. 
Such a situation cannot be exactly duplicated, though, due to the nature of acoustic waves 
to spread and travel in all directions. In order to restrict the acoustic signal to only one 
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direction, it is proposed that the sound field be set up in a cylindrical, resonant sound 
chamber so that the analysis is properly limited to axial wave modes only. 

B. FUNDAMENTAL mEAS 

Unlike the well known dependence on the Reynolds and Prandtl numbers found in 
conventional cross flow over bluff bodies, the issue of heat transfer in the presence of an 
acoustic field is significantly more complicated due to the presence of a multitude of 
competing length (or time, or velocity) scales. The ways in which these length scales may 
be ordered are many and lead to numerous distinct parameter regimes with quite drastically 
different flow properties, and hence heat transfer properties. In order to closely examine the 
properties of heat transfer in oscillatory flows then, it is first necessary to enumerate some 
of the different parameters and variables involved. By establishing and maintaining a set of 
criteria surrounding these parameters, and by modeling the test apparatus to conform to them, 
the job of evaluating different regimes of flow and drawing correlations from the data 
obtained will become possible. 

This section lists the most important of these parameters and provides reasoning for 
the choices involved, keeping in mind the desire to maintain the most basic and core set of 
conditions that flow around the test cylinder be laminar, incompressible and attached, and 
that only the steady-state solution be considered. 

1. Criterion A 

Following the first assumption that the acoustic streaming flow around the test 
cylinder be incompressible, two different length scales must be specified. The first of these 
requires that the relative size of the test cylinder be very small compared to the radian 
wavelength of the acoustic field, where the radian wavelength is defined as 


and the characteristic length scale of the test cylinder is chosen as the radius, a. The radian 
wavelength can be related to the frequency by: 
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( 2 ) 


X - — - s!L - ^ 

2 % 2 % (0 

Now, by asserting thata«2, and designating the ratio between the two as it follows that; 



(3) 


and finally that 



(4) 


The criterion X ^ ensures that radiation effects due to the acoustic streaming are 
negligible, as presented by Lighthill (1963) and indicates that it is only the local acoustic 
field conditions that are of importance. 

2. Criterion B 

The second criterion which is required to support the assumption that flow be 

incompressible is derived from the relationship between displacement amplitude of particle 

— ~I 

oscillation in the sound field, A , and the cylinder radius, a. The ratio between the two, —, 

a 

will be designated as the amplitude parameter, e, and dictates whether or not separation will 
occur. 

When the amplitude parameter is very small, 

€ = - « 1 (5) 

a 


the particles in the sound field move a very short distance along the cylinder wall before 
reversing their direction. This ensures that the flow remains attached, with little chance of 
separation occurring, and hence the flow will remain laminar at all times. It can also be 
noted that £is directly proportional to the pressure ratio Pq / and can therefore take on 
much larger values for strong acoustic fields. This is accomplished by observing that the 
displacement amplitude of a particle in the flow is directly related to its velocity and that the 
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particle velocity is in turn related to the pressure ratio by the following (in a plane standing 
sound field) 



b) 


and 


^0 = 


cPq 


( 6 ) 


(7) 


therefore 


—1 


( 8 ) 


Yet another form of this parameter often used in the literature is called the Keulegan 
Carpenter number and is defined as KC=UJ2af and can additionally be expressed as 


KC = ne (9) 

Of importance is the fact that the product of the parameters defined in the first two criteria 
(A and B) is the flow Mach number, which can be defined as 

^0 

M=xe=— (10) 

c 

When criteria A and B are satisfied, M«l, and this in turn is the second condition which 
satisfies the assumption of incompressible flow. 

3. Criterion C 

The Stokes boundary layer thickness S is related to the kinematic viscosity and the 
radian frequency of oscillations by 
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and is the well known length scale which is a measure of the extent of viscous effects in an 
oscillatory flow. A frequency parameter can be defined as follows 


A 2 = («)2 ^ 

6 V 


( 12 ) 


For the case when A^»l, the Stokes shear layer is confined to a narrow region and the 
acoustic streaming effect appears as a slip velocity along the cylinder surface. Utilizing the 
knowledge that the boundary layer thickness is on the order of 103 and imposing the 
condition (somewhat arbitrarily) that 


it follows that 


106 


(13) 


> 1600 


(14) 


is a good criterion to ensure “large” values of the frequency parameter. 

The frequency parameter may also be often found in the literature in the form of 
P =(2aff/\, and can be expressed as 


P = (-)A2 (15) 

4. Criterion D 

When criteria A - C are satisfied, the acoustic streaming velocity is of magnitude 
O(eUo). A streaming Reynolds number, i?,., can then be defined as 



( 16 ) 
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Through substitution for e and Uq 


which yields 




aci V 



cov 


—(—)' 
yPm 


Rs = e'A' 


(17) 


(18) 


This streaming Reynolds number becomes the driving factor in determining what will be the 
primary mode of heat transport within the region due to the acoustic streaming flow. Stuart 
(1966) demonstrated that when R/<1, a Stokes flow becomes prevalent in the outer region 
while a boundary layer flow is predominant when this parameter takes on values much 
greater than one. In order to ensure that forced convective heat transfer is dominant then, we 
impose the following constraint 

R, » 1 (19) 


5. Criterion E 

It was first observed by Honji (1981) that flow around a cylinder will become 
centrifugally unstable and separate into vortices as the amplitude of particle oscillation 
increases. This instability occurs in the Stokes layer where the flow is parallel to the direction 
of particle oscillation. This was confirmed by Hall (1984) who conducted a linear stability 
analysis on the unsteady boundary layer in the high-frequency limit. He found that a critical 
value of the Reynolds streaming number exists for which instabilities begin to form, namely 
when jRj becomes greater than 4.24A. Recent work in the area by Sarpkaya (1986) provides 
further explanation for this phenomenon. Since vortex shedding can make a large impact on 
the heat transport from the cylinder, we will maintain the following criterion, which can be 
expressed in two different ways as 


R^ < 4.24A 


( 20 ) 
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or, alternately as 


v/A 

An example of what these vortices may look like is shown in Figure 6. 

6. Criterion F 

In order to maintain the condition in which there is minimal influence from the 
buoyancy effects of natural convection as compared to the forced convective heat transfer 
due to the acoustic streaming effects, a dimensional analysis of the governing equations 
produces the following requirement for the Grashof number 



( 22 ) 


For the case when Gr / ~ 1 or greater, buoyancy effects must be taken into account and 

any heat transfer correlations developed will have to be modified accordingly. 



Figure 6. Vortices Due to Oscillatory Flow Around a 
Cylinder. (Sarpkaya, 1986) 





7. Criterion G 

This final criterion stems from the desire to maintain only a plane standing sound 
field within the test chamber. This means that only the axial wave components excited by 
the resonant acoustic signal be present, without interference from transverse modes such as 
the azimuthal modes and radial modes. It may be shown from theory that there is a certain 
cut-off frequency developed from a solution of the wave equation, below which the 
transverse modes will not be present. For the cylindrical waveguide geometry being used, 
the solution for the above condition is obtained in the form of appropriate roots of Bessel 
functions. If L is the length of the test chamber, the longitudinal frequency can be expressed 
as 



(23) 


where I is the F‘, 2"'^,..., mode number. The transverse frequency mode is expressed by 


j: ' 

f = (X 

Jmn mi 


tzD 


(24) 


where D is the test chamber diameter and represents the eigenvalues obtained from roots 
of the Bessel functions. Transverse modes will be present if /^ = , so it is desired to 

maintain/; well below^ and/ now represents the maximum frequency possible while still 
maintaining this criterion. By substituting for both frequencies, the condition becomes 


lc_ 

2L 


< 


cc„ 


c 

%D 


which is reduced to 


/ 




a„ 


% 


By introducing a new parameter called the aspect ratio. 


(25) 


( 26 ) 
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( 27 ) 



and finding the smallest root of this Bessel function, 

/ 

—^min = 0.586 
n 


(28) 


the criterion now becomes 


— < 0.586 
IZ 


(29) 


which can be rearranged to finally get 


where 


Z = 


D 


> 0.85 L 


(30) 




f 


(31) 


This gives a relationship between the geometry and the maximum frequency. But since the 
maximum frequency is already defined using criteria A - F, criterion G in fact gives us the 
maximum diameter that can be used, or if the diameter is also given, it determines the length 
of the chamber instead. 

8. Basic Flow Description 

When a cylinder is immersed in a standing acoustic field and all of the previous 
criteria have been met, particle oscillations will initiate the most basic form of an acoustic 
streaming flow. This steady flow (Figure 7) is symmetrical about two axis, is circular in 
nature, and includes a well defined boundary layer. The boundary layer is actually quite small 
and is greatly exaggerated in the figure for clarity. Although this is not the only flow which 
is present, it does represent the flow pattern that has the biggest impact as a heat transport 
mechanism when these criteria are satisfied. 
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Figure 7. Outer Acoustic Streaming Flow (also known as Boundary Layer Flow) 


As discussed earlier, this fluid flow regime is the main focus of this investigation in 
which at first we consider only laminar, unseparated flow around the test cylinder. One of 
the goals of this experimental study is to determine the range of values for the previously 
discussed parameters where one flow regime transitions to another and relate them in terms 
of the streaming Reynolds number. It is believed that such transitions will be reflected in the 
heat transfer behavior. 

In order to define a specific range for laminar flow, though, these criteria needed to 
refinement by asking the questions; how small is “much less than one” or how big is “much 
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greater than one”? It was determined that Criterion A in Eq. 4 would be met by picking 
X < 0.1 and Criterion B in Eq. 5 would be met by choosing e <0.3. This ensures 
incompressible flow around the test cylinder and that flow remains attached. Criterion C 
becomes > 1600, confining the Stokes shear layer to a narrow region. The above ranges 
have been fixed based on prior experience and preliminary experimentation, but are not by 
any means intended to be “hard and fast”. They may indeed be modified if necessary. 
Criteria D and E require verification after each experimental run. 

C. APPARATUS 

As previously discussed, the experiment consists of a heated test cylinder placed 
normal to a simple acoustic standing wave within a resonating sound chamber (Figure 8). 
This general description can be broken down into three major components; the test cylinder, 
the sound chamber and the acoustic electronics package, which are further described below. 



Figure 8. Experimental Test Apparatus. 
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1. Test Cylinder 

It was originally proposed that a Watlow stainless steel cartridge heater be used as 
the test cylinder for the experiments. They are available in varying diameters and 
additionally feature an imbedded type “J” thermocouple placed at the midpoint of its length. 
This presented problems, however, since the heating along the length of the cartridge heater 
was uneven and the relatively large thermal resistance of the stainless steel produced large 
variances in surface temperature along the length. Since constant surface temperature is a 
feature which is very important to the analysis of heat transfer characteristics, this was 
unacceptable. Therefore, a copper sheath was designed to fit over the cartridge heater. It 
was assumed that the large thermal conductivity of copper would even out the axial surface 
temperature gradient. In addition, silicon oil was liberally applied to the inside of the copper 
sheath prior to insertion of the cartridge heater before each experimental run (Figure 9). Ihis 
provides better thermal contact in the narrow annular gap (approximately 2 mm ) between the 
sheath and heater, preventing air gaps which can cause local temperature discontinuities at 
the surface. This arrangement (as was later verified. Appendix B) provides for an axial 
temperature variation of less than 0.5°C from one end of the heater to the other for the 
temperature range of interest in this experiment. A picture of the test cylinder is shown in 
Figure 10. 
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Figure 10. Photo of Test Cylinder inserted into the Sound Chamber. 








As was stated earlier, the surface temperature of the test cylinder is a required datum 
point for the analysis. Since the only provision for temperature measurement is at the 
thermocouple placed at the center of the test cylinder, the equivalent resistance of the steel/ 
oil/copper circuit is required in order to deduce the surface temperature from the cartridge 
heater center temperature as measured by the imbedded thermocouple. Appendix B gives 
a detailed analysis of the derivation of the resistance, which is approximately 1.019 K/W. 

The cartridge heater receives its power from a Kikusui Model PAR 160A regulated 
DC power supply. It provides power control measurement down to 0.01 amps and 0.01 
volts. Thermocouple measurements are provided by a Keithley Model 740 scanning 
thermometer system. Calibration data for all thermocouples and the thermocouple reader is 
provided in Appendix B. 

2. Sound Chamber 

The purpose of the sound chamber is to provide an environment through which 
acoustic signals of various frequencies can be used to excite a resonant standing wave. In 
order to accomplish this for all frequencies which may be used during the experiment, a 
resonant chamber that would be adjustable in length was highly desired. Figure 11 shows 
the final configuration of the test chamber while Figure 12 shows a photo of the test 
apparatus. 



Figure 11. Sound Chamber Assembly. 
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Figure 12. Photo of Test Apparatus. Acoustic Driver is on to the right. 





The test chamber itself is a plexiglass tube approximately two meters long which is 
mounted firmly to a plexiglass base plate. The acoustic driver is mounted at one end and 
provides the source of the acoustic signal. The opposite end has a movable flat end plate 
which has an o-ring seal to help isolate the acoustic wave in the chamber. The end plate 
serves as the rigid end termination while the o-ring provides good sound confinement, as 
well as good stability for the end face so that it doesn’t become offset to either side as its 
position is varied along the length of the enclosing resonant chamber. 

A hole through this end plate provides access for a pressure transducer which is used 
to help deduce the pressure ratio, and hence the sound pressure level for that configuration 
of frequency and input power. The pressure transducer used is an Endevco model 851 OB-5 
which has an output in millivolts and has a pressure sensitivity of 50.89 mV/psi. This is 
connected to a preamplifier with a gain of 100 and provides output to both an oscilloscope 
and a Hewlett-Packard model 34401A multimeter. Measurement of the voltage output from 
the multimeter is important to determine when resonance has occurred within the chamber 
for as the frequency is varied, the output voltage from the microphone decreases on either 
side of the resonant operating point. 

The oscilloscope provides a visual representation of the time trace of the acoustic 
signal at the end of the sound chamber, corresponding to a pressure antinode, and is used to 
ensure that the signal remains sinusoidal throughout the experimental range of powers and 
frequencies. During the initial stages of the experiment, the oscilloscope allowed for the 
discovery of interference patterns in the sinusoidal waveform as caused by higher order 
harmonics at high SPLs (> 155 dB). In order to limit the interference that was present, the 
use of “sleeves” within the chamber was recommended to detune the resonant mechanism. 
This would prevent the harmonics from being integral multiples of each other and thereby 
prevent them from reinforcing each other to form interference patterns. By placing a sleeve 
at an appropriate spot in the chamber, some of the high frequency harmonics leading to 
interference could be eliminated, allowing for even higher SPLs to be achieved before 
interference occurred. 
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3. Acoustic Electronics Package 

Of great importance to the experiment is the ability to generate a nearly pure 
sinusoidal waveform at varying frequencies and high power ranges. As the strength of the 
signal generated increases, the effect of the flow around the test cylinder becomes more 
pronounced, enhancing the heat transfer characteristics of the system. 

The acoustic signal being generated at the driver end of the sound chamber is 
provided by a Hewlett-Packard model 33120A arbitrary waveform generator. It sends a 
sinusoidal waveform at the proper frequency through a Techron model 7540 power supply 
amplifier to a JBL model 2490H acoustic compression driver. 

A more detailed review of each item of equipment used in the system is provided in 
Appendix D. 

D. EXPERIMENTAL METHOD 

Prior to gathering data for this experiment, it was first necessary to develop a 
coherent plan with which to approach the problem. The first step required was to find 
specific frequencies at which resonance occurred within the chamber, and which would also 
provide a velocity anti-node at the position of the test cylinder. By looking at the geometry 
of the problem (Figure 13), it became obvious that these limiting factors could be met by a 
combination of adjusting the end plate distance from the test cylinder, as well as the signal 
frequency, so that the length L from the heated cylinder to the end plate termination was an 
odd multiple of XIA. This could be further refined to state that the value of 4Lf/c needs to be 
an odd integer. When this condition was satisfied, the requirement of a velocity anti-node 
occurring at the cylinder location was met. 

The next thing needed was an estimate of the maximum pressure ratio, and therefore 
the sound pressure level that could be obtained. This was achieved by increasing the 
amplitude of the input signal waveform until the output waveform on the oscilloscope began 
showing traces of interference or other disturbances. When the maximum input amplitude 
was obtained, the multimeter output was recorded. 
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where n = 5 
for this case 

X — 



Figure 13. Proper Geometry so that Maximum Velocity occurs at the Test 
Cylinder. 


In order to convert the output voltage to an actual pressure ratio and sound pressure 
level, it was necessary to first understand what form the multimeter uses to present the 
output. The multimeter gives the voltage output in terms of the true RMS value of the 
sinusoidal signal, or rather VJsjl, where Vo represents the voltage amplitude from zero to 
peak of the sinusoidal signal. Recall that the true RMS value for the pressure can also be 
expressed as P^\f2. Then, the following relationship holds true, that 

(32) 

where S is the sensitivity of the pressure transducer in mV/Pa. 

The sensitivity of the pressure transducer after passing through the 100 gain setting 
of the preamplifier was converted and expressed as 0.74 mV/Pa. By substituting into the 
above equation, the pressure amplitude becomes 
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p - 

® 0.523/nWPa 

where it is noted again that ! \fl is the multimeter output. 
From this, a pressure ratio can be defined as 


( 33 ) 


PR = 


(34) 


where is the mean ambient pressure of 101 kPa. The pressure ratio is typically expressed 
in terms of a percentage. 

The sound pressure level is defined as the logarithmically scaled ratio of the RMS 
pressure and a predetermined reference pressure, P^f, where Pref is chosen to be 20 pPa for 
gases (by convention). This gives 


SPL = 20 


log 




10 


ref 


(35) 


In order to obtain an idea of how strong the acoustic signal is during the experiments. Table 
1 gives the sound pressure level for various activities 


Activity 

Sound Pressure Level 

Pressure Ratio 


(dB) 

(%) 

Normal Conversation 

60 

<0.001 

Jet Airplane at Take-off 

90 

<0.01 

Pain Threshold 

120 

0.28 

Minimum Experimental Level 

150.6 

0.94 

Maximum Experimental Level 

161.2 

3.20 


Table 1. Samples of SPL and PR for Comparison. 


After finding the pressure ratio, and therefore the SPL, the actual data extraction 
phase of the experiment could be initiated. It was preceded by with an understanding of what 


29 



















type of data was needed for analysis. A crucial element of this experiment is the need to 
obtain a broad base of data with which to incorporate the results. In order to obtain this, 
experiments were run at several different pressure ratios for each resonant frequency found, 
starting from as low as 0.9% and building up to the maximum pressure ratio by increments 
of 0.1%. In order to obtain a spread of data points at each pressure ratio evaluated, the test 
cylinder was heated to approximately 8, 12 and 16 degrees above the ambient temperature. 
This kept the power requirements low and reduced the amount of thermal input into the 
ambient air within the sound chamber. This latter point was necessary as the ambient 
temperature could rise as much as 0.3 degrees during a single experimental trial. In order to 
ensure reproducibility and reduce the effect of anomalistic behavior, three runs at each 
temperature point were conducted to ensure consistency of data. The selected temperatures 
were only guidelines and were not meant to be hard set points for the experiment. Instead, 
they were treated as aim points with an acceptable range of ± 1 degree. Therefore, in order 
to produce an even broader spread of data, the power input to the test cylinder was varied 
slightly for each of the three trials at each specific temperature point. 

Once the selected frequencies and pressure ratios had been determined and a suitable 
starting pressure ratio and temperature had been obtained, the experimental process was 
initiated. In order to obtain the selected pressure ratio for a particular set of runs, it was 
necessary to get an idea of the settings required for each specific piece of equipment. This 
was done by selecting the appropriate frequency on the waveform generator and modifying 
the power amplification on both the waveform generator and the power amplifier until the 
appropriate pressure ratio was obtained from the multimeter. The frequency was then 
adjusted to fine tune the resonance. A check of the oscilloscope at this point ensured that the 
signal being generated was of the right waveform and that interference was not occurring. 
Then the power amplifier to the acoustic driver was turned down to zero after noting the 
level at which it was set. This allowed for obtaining the correct pressure ratio in a quick 
manner by simply turning the power amplifier up to the previously noted value. 

It was necessary at this point to ensure that the test cylinder was properly prepared. 
This entailed introducing approximately ten drops of silicon oil into the copper sheath and 
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inserting the cartridge heater. The cartridge heater would be completely immersed in the 
silicon oil after being fully inserted. After a period of time, there would be some loss of 
silicon oil due to the wick action of the thermocouple and power leads emerging from the 
top. This was insignificant during the mns required for a single pressure ratio and caused a 
negligible change in the center-to-surface resistance as noted in Appendix B, but it became 
good practice to add some of drops of oil each time a new set of mns at a different pressure 
ratio were to be taken. 

The test cylinder was then inserted into the sound chamber, making sure that the 
bottom of the cylinder was resting in a shallow indentation specifically machined into the 
inside face of the chamber. Power to the cartridge heater was then turned on and set so that 
a steady state temperature of approximately eight degrees above ambient would occur when 
the acoustic signal was present. This became more of an art form and required familiarity 
with the system to accomplish it with any degree of accuracy. 

Since the power to the acoustic driver at this point was still at zero, the power to the 
cartridge heater would drive the temperature of the test cylinder past the projected steady 
state temperature. As it approached the projected temperature, though, the power to the 
acoustic driver was increased to the previously noted set point. This provided the fluid flow 
at the predetermined pressure ratio, in effect beginning to cool the test cylinder until it 
reached a steady state condition. At this point, the rate of energy input to the cylinder was 
equal to the rate of energy being convected away from the cylinder. By monitoring the 
interior thermocouple temperature, it was easy to see when the lowest temperature was 
reached. When the temperature began to rise once again, it was determined that the steady 
state condition had been reached and that the resultant temperature rise being witnessed was 
due only to the test cylinder transferring heat into the surrounding fluid medium, thus raising 
the overall ambient temperature. 

Two situations other than the ideal one presented above occurred frequently due to 
the coarse means of trying to arrive at the desired cylinder temperature. If upon engaging the 
power to the acoustic driver the temperature of the test cylinder did not decrease but 
continued to increase at a slower rate instead, then the power being supplied to the cartridge 
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heater was deemed too high and the voltage reduced until a decrease in temperature was 
witnessed. If upon engaging power to the acoustic driver the temperature of the test cylinder 
were to continue decreasing lower than the desired temperature range, it was an easy matter 
to increase the voltage to the cartridge heater to increase the steady state temperature 
solution. 

When the steady state solution point was reached, the frequency and microphone 
voltage were noted, as was the voltage and current supplied to the cartridge heater and the 
temperature of the thermocouple in the cartridge heater. The power to the heater and the 
power to the acoustic driver were then simultaneously turned off and the test cylinder was 
quickly removed from the sound chamber. A second thermocouple was then introduced 
through the access hole in the sound chamber (from which the test cylinder had just been 
removed) such that the location of the thermocouple was approximately at the center (i.e., 
along the axis) of the sound chamber. This thermocouple temperature was then monitored 
until it “plateaued” and provided the measurement of the ambient temperature within the 
sound chamber. 

This completed a single experimental trial for a specified resonant frequency, 
pressure ratio and temperature. The procedure was carried out a total of nine times for each 
different pressure ratio. Once the six outputs for each run were recorded, they were then 
transferred over to a spreadsheet where all other significant parameters were computed 
automatically. The following section delineates the various calculations performed in the 
spreadsheet. 

E. EXPERIMENTAL CALCULATIONS 

The first calculation desired from the spreadsheet entails finding the actual power 
being supplied to the cartridge heater. This is expressed in terms of the current and voltage 
outputs from each experimental ran. 


P = IV 


(36) 
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Once the power to the cartridge heater was known, the surface temperature of the test 
cylinder itself could be calculated. Knowing the resistance of the thermal circuit as given in 
Appendix B, the difference between the interior temperature and the surface temperature is 
simply 


Ar = (37) 

Utilizing the interior temperature output as provided by the thermocouple embedded in 
the cartridge heater, the surface temperature of the cylinder is then expressed as 


T. - PR 

c eq 


(38) 


Once the surface temperature is known, the difference between it and the ambient 
temperature, as measured during the experiment, is calculated. The result is then combined 
with the power calculated in Eq. 35, as well as with the external surface area of the test 
cylinder, to find the convective heat transfer coefficient 


h = 


P 

MT, - TJ 


The Nusselt number is then derived from the following equation 


(39) 


Nu = — (40) 

k 

where d is the test cylinder diameter and k is the thermal conductivity of air. 

The next step is to calculate the various criteria as previously listed in the theory 
section. In order to find the length scale ratio, %, the speed of sound within the chamber must 
first be calculated using the ambient temperature 


c = ^yR(T^ H- 273.15) 


(41) 
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The value of % can now be found from Eqs. 1 - 4. The amplitude parameter, e, is derived in 
the spreadsheet by combining the pressure ratio from Eq. 33 along with Eq. 8. The third 
parameter calculated is the frequency parameter, A^, from Eqs. 10 and 11. 

The streaming Reynolds number, R^, is found by using Eq. 16. This, however, does 
not represent the true value at the test cylinder position due to its not being precisely at the 
velocity antinode. A corrected value for R^ can be deduced by finding the particle velocity 
offset between the velocity antinode location and the test cylinder position. This offset is 
derived by knowing the frequency of the sinusoidal signal being generated, the speed of 
sound within the chamber from Eq. 40, and the distance from the test cylinder to the end 
plate termination face, L. Hence 


17o|sin( — - )| 

c/f 2 ' 


or 



R, = /fjsin(^)|^ 

‘^corrected ^ ^ 


{‘m which it was found that the corrected value for R^ had less than a 0.1% error due to the 
slightly displaced location of the cylinder. 
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IV. RESULTS AND DISCUSSION 


Nearly 600 experimental trials were performed throughout the course of this study 
which produced results for 183 distinct data points. Data were obtained for five different 
frequencies at various pressure ratios ranging from 0.9 % to 3.2 %. Three separate series of 
trials were performed at each acoustic signal setting, (i.e., at each frequency and pressure 
ratio setting) corresponding to three separate settings for the driving temperature difference 
between the test cylinder surface and ambient conditions. The resultant values of the 
streaming Reynolds number ranged from 40 to 1070 while the corresponding Nusselt 
numbers obtained varied from 8 to 38. Figure 14 is a parameter map as suggested by 
Richardson (1967) which shows the range of values for the experimental data covered plotted 
as a function of the amplitude parameter versus the frequency parameter, and delineates the 
expected different regimes of flow. The data obtained cover a very narrow regime of this 
parameter map as was intentionally planned for this experiment. Now it can be seen that the 
heat transfer results obtained through experimentation are quite evenly distributed between 
two distinct regions on the map. Region A represents the regime in which the flow is 
expected to remain laminar, incompressible and attached, and outer acoustic streaming is the 
main heat transport mechanism. This region is well understood in theory but has yet to be 
thoroughly verified through experimentation. It is anticipated that the heat transport 
mechanism for the data in region B will be presented as a combination of effects, including 
that of vortex shedding, as predicted by Honji (1981), Hall (1984) and Sarpkaya (1986). 

The heat transfer results are presented as plots of Nu vs. Rj (Figure 15). Here the 
difference between the two regions of varying heat transport mechanisms becomes more 
apparent. A break clearly occurs in the region where R^ ~ 240 (± 20) and it can be concluded 
that there is some critical point in this range where there is a transition in the flow at which 
vortex shedding begins to become a dominant factor in the heat transport away from the test 
cylinder. In order to better examine the differences between these two regimes, the data is 
divided into their respective groups and individually analyzed. 
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Figure 14. Parameter Map of Expected Heat Transfer Regimes as presented by 
Richardson (1963): Convection by Inner Acoustic Streaming (A), by Outer Acoustic 
Streaming (B). 
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Figure 15. All Data Plotted as a Function of Nusselt Number vs. Streaming Reynolds 
Number. 
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A. LAMINAR, ATTACHED FLOW REGIME 

Figure 16 is a plot of the heat transfer results in terms of the Nusselt number versus 
the streaming Reynolds number for the data in which < 240, and for which criteria A 
through C (Eqs. 4,5,15) have been met. Those data points which do not meet these criteria 
have been disregarded. The critical parameters for the remaining data are as follows: 

X<0.1 

e < 0.3 

A2>1800 

R 

< 4.5 

A 

Theory clearly indicates that the dependency of the Nusselt number on the streaming 
Reynolds number in this regime is of the form Nu = xPr^s®’^. Since the Prandtl number 
remains constant throughout the experiment, the solution for this dependency be can further 
simplified as Nu = CR^” ^ where the term “C” encompasses both the Prandtl number and the 
qualitative constant of the previous equation. However, this solution form is only valid for 
“large” values of R^. Since the theory does not provide a definite limit for what qualifies as 
“large”, this criterion had to be determined from a careful examination of the experimental 
data. From Figure 16, it can be observed that there is indeed a break point at R^ ~ 130 where 
the results diverge into two separate solutions. It was found that the square-root dependency 
on Rj does not significantly change past this value, and it is therefore suggested that this may 
be in the range of the lower limit of “large” values of Rj. For those values in which Rs> 130, 
a curve fit of the heat transfer characteristics results in a solution of the form 

Nu = 0.94R,®^ (45) 

and it can therefore be determined that the range of values 130 < R^ < 240 is representative 
of “large” values of the streaming Reynolds number. The values below R^ = 130 are 
excluded as not being large enough due to a variety of reasons as described later. 
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Davidson (1973), in an extension of the work by Richardson, analytically and 
numerically tackled the problem of heat transfer from a cylinder in a strong acoustic field in 
great detail. He obtained a correlation of the dependency of the Nusselt number on both the 
Prandtl and streaming Reynolds numbers. The correlation, as extracted from the work of 
Davidson by Gopinath and Mills (1993) for this regime, is of the form 

Nu = 1.388PrO-^2R/-5 (46) 

By taking Pr = 0.7 for air for these experiments, the equation then becomes 

Nu = 1.07R/-^ (47) 

The experimental fit in Eq. 45 under predicts by about 13%, but supports this correlation 
well within the limits of uncertainty. 

Figure 17 is a plot of the region in which Eq. 45 is valid, and includes the 
experimental uncertainty of each data point as derived in Appendix C. It can be observed 
from this plot that the deviation from the curve fit in this range of values for R^ is well within 
experimental uncertainty limits. 

Although the lower end of the range for “large” R^ in which the predicted solution is 
valid has been determined to be 130 for the results obtained from these experiments, it is by 
no means an absolute boundary. Even though the resultant heat transport characteristics 
deviate significantly below that point in the region where “intermediate” values of Rj are 
present, there are several factors which could account for part of the discrepancy, especially 
in the region around 100 < R^ < 130. These include uncertainty due to equipment limitations 
and the effects of natural convection and conduction of heat away from the test cylinder. 

The effect of natural convection on most of the intermediate values of is 
negligible, though, as characterized by the very low ratio of the Grashof number to the square 
of the streaming Reynolds number (except at very low values of Rj). The effect due to 
conduction is much harder to quantify in so simple a form, though. Conduction can and does 
occur during the experiment at two separate places where the test cylinder is in contact with 
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other apparatus. The top of the test cylinder is in contact with the “plug” used to hold it in 
place, while the bottom of the test cylinder is allowed to rest in a small depression on the 
wall of the plexiglass sound chamber. Realizing, though, that the contact area at both points 
is very small, and that the thermal resistance of the plug and the plexiglass wall are both 
relatively high (due to their low thermal conductivity), the heat transport away from the test 
cylinder can be assumed to be negligible as well. 

The uncertainty in the results due to the equipment limitations, though, does have a 
very significant impact on the results. The calculated value of the Nusselt number has an 
uncertainty of up to 20% (see Appendix C) depending upon the power being dissipated by 
the test cylinder. Since the experiment revolves around finding the heat transfer 
characteristics at specific values of temperature difference between the test cylinder surface 
and the ambient conditions, correspondingly low power dissipation from the cylinder occurs 
as the streaming Reynolds number decreases. Therefore, the region of intermediate values 
of Rj have relatively low electrical heat dissipation, and hence low current values associated 
with them, dropping to as low as 0.06 amps in some cases. Since the equipment uncertainty 
for the current reading is of the order of the last digit present, this particular component of 
the Nusselt number has an uncertainty of nearly 18% by itself and greatly influences the 
overall uncertainty. Therefore, it may be more accurate to define the lower limit of large R^ 
values as somewhere between 100 and 130. However, the error due to equipment 
measurements is not enough to compensate for the disparity between theory and experiment 
at values much less than 100. 

B. SEPARATED FLOW REGIME 

The second regime which this experiment encompasses is that in which vortex 
shedding and other forms of unsteady flow begin to affect the heat transport characteristics. 
Figure 18 is a plot of the resultant data in this regime as obtained during experimentation in 
terms of the Nusselt number versus the streaming Reynolds number. A curve fit of the data 
results in a solution of the form 


Nu = 0.31R,°-*^ 


(48) 
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Figure 18. Unstable Regime. 
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There is no theory for this flow regime with which to compare the results, but it is reasonable 
to assume that this curve fit is representative of the correct solution. As expected, the 
unstable flow resulting at higher values of would increase the heat transfer rate from the 
cylinder, and hence the stronger dependency that the Nusselt number has on the streaming 
Reynolds number, as opposed to the usual square root dependence. 


44 



V. CONCLUSIONS 


Experiments were conducted to observe the convective heat transfer rates to an 
isolated cylinder in an acoustic standing wave. A comparison of the various length scales 
and other parameters was conducted, and the experimental method stated. During the 
experiment, the properties of the acoustic field were varied to provide a large base of data 
which was then analyzed and discussed. Several regimes of interest were investigated and 
the results presented. Figure 19 is a plot of all data obtained with curve fits for the regimes 
of interest. 

Essentially, heat transfer from a cylinder in a zero-mean oscillatory flow as 
represented by an acoustic standing wave can be divided into at four separate regimes in 
which different heat transport mechanisms dominate. For very low values of the streaming 
Reynolds number, R^, convective effects due to the acoustic field are negligible and natural 
convection is then the dominant mode of heat transport. For intermediate values of Rj, there 
is a stronger dependence on R^ but not yet on the order of R^ since Rj is still not large 
enough for flow to be of the boundary layer type. Buoyancy effects are comparable in the 
lower end of this regime, becoming small for larger values of R^. For much larger values of 
Rj, past 100 or so, an acoustic streaming flow presents itself in the boundary layer, resulting 
in a square-root dependency on Rj. Experimentally, the results obtained in this regime 
closely match the expected theory, and the heat transfer characteristics may be estimated by 
Eq. 45. Finally, past a critical value of R^ -240 (which confirms well with theory), an 
unstable flow with vortex shedding begins to take place at the surface of the cylinder, 
increasing the dependency on R^ which the overall heat transfer solution has. The heat 
transfer characteristics in this regime may be estimated from Eq. 48. It is these last two of 
the above regimes that formed the focus of this study. 
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VI. RECOMMENDATIONS 


Several directions for further study may be suggested by the data obtained from this 
work. One factor which may need further investigation is a study of the effects of varying 
the aspect ratio of the chamber diameter to the cylinder length to determine whether it was 
large enough for these experiments. A ratio of 15 was used for this experiment, and was 
assumed to be large enough to discount flow effects caused by the walls of the chamber, but 
this number was chosen somewhat arbitrarily and only a detailed experimental study can 
determine the actual effect. 

Another possible source of error which cannot be accurately accounted for concerns 
itself with the experimental method used. Specifically, the measurement of the ambient 
temperature within the chamber at the time of each experimental trial. The method used, that 
of removing the test cylinder from the sound chamber after its temperature has been recorded 
and the acoustic signal has ceased, then placing a thermocouple through the hole it has just 
evacuated in the chamber wall, can be improved upon. During the time required while 
waiting for the thermocouple temperature to peak, the heated air in the sound chamber is 
rising to the top of the chamber and the accuracy of correlating the thermocouple temperature 
reading to that of the ambient temperature at the time of the experimental trial is left in 
doubt. A better method would be to permanently affix the thermocouple in the chamber so 
that simultaneous measurement of both the cylinder temperature and the ambient temperature 
can be taken. 

Additional research into three different areas of the problem come to mind. First, a 
test of the effects of placing the test cylinder horizontally in the sound chamber is suggested, 
although there should be little if any dependency on this orientation since natural convection 
effects are small for the strong acoustic fields being used. Completely new geometries may 
also be tested which would mimic actual heat exchanger component shapes expected in a 
thermoacoustic engine. Finally, additional research using different gases in the sound 
chamber would provide data on the dependency of the Nusselt number on the Prandtl 
number. 
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APPENDIX A. CALIBRATIONS AND CALCULATIONS 


Several pieces of the experimental apparatus required some form of calibration, or 
calculation of a specific parameter, prior to initiating the experiments. The most important 
equipment items of concern were the “unattached” J-type thermocouples which were to be 
used for various applications throughout the experiment, as well as the thermocouple which 
was embedded in the cartridge heater. In addition, an equivalent thermal resistance for the 
cartridge heater/silicon oil/copper sheath circuit needed to be calculated along with 
assurances that the linear temperature distribution along the test cylinder was within 
reasonable limits. An additional study was performed to analyze the heat transfer effects on 
the test cylinder due only to natural convection. 

Three J-type thermocouples were used throughout the experiment and were the first 
items to be calibrated. Since accurate temperature information was crucial to the reliability 
of the data obtained through experimentation, the thermocouples were tested to see if any of 
them showed a tendency to read either higher or lower than the actual temperature (a 
somewhat common occurrence). The embedded thermocouple in the cartridge heater was 
also tested for the same reason. All thermocouple leads were attached to the thermocouple 
reader to be used throughout the experiment to ensure that the entire circuit was tested 
concurrently. 

The unattached thermocouples and the cartridge heater were all placed in an ethyl 
glycol solution belonging to a Rosemont Model 913A calibration bath. A Rosemont Model 
920A commutating bridge, which utilizes a precision temperature probe as its input, provides 
the reference temperature. The temperature of the bath was then set at varying points 
between 22° and 48°C, the expected experimental temperatures being well within that range. 
Table 2 shows the results of the calibration. All thermocouples were found to read within 
0.1 °C of the reference temperature for all cases. This was well within the possible 
uncertainty of 0.5°C listed for J-type thermocouples. 

Next, Figure 20 shows the thermocouple arrangement which was used to 
experimentally derive the equivalent thermal resistance and the linear temperature 
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Reference 

Cartridge 

Thermocouple 

Thermocouple 

Thermocouple 

Maximum 

Temperature 

Heater 

“A” 

“B” 

“C” 

Deviation 

(C) 

(C) 

(C) 

(C) 

(C) 

(C) 

22.83 

22.8 

22.8 

22.8 

22.8 

<0.1 

23.96 

23.9 

23.9 

23.9 

23.9 

<0.1 

26.83 

26.8 

26.8 

26.8 

26.8 

<0.1 

29.63 

29.6 

29.6 

29.6 

29.6 

<0.1 

32.49 

32.5 

32.5 

32.5 

32.5 

<0.1 

35.37 

35.4 

35.4 

35.4 

35.4 

<0.1 

38.04 

38.1 

38.0 

38.0 

38.0 

<0.1 

40.74 

40.8 

40.8 

40.8 

40.8 

<0.1 

43.02 

43.0 

43.0 

43.0 

43.0 

<0.1 

45.68 

45.7 

45.7 

45.7 

45.7 

<0.1 

48.19 

48.2 

48.2 

48.2 

48.2 

<0.1 


Table 2: Thermocouple Calibration Data (all values in °C) 



Figure 20. Thermocouple Arrangement for Calibration Tests. 
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distribution data. The three unattached thermocouples which were previously calibrated 
were securely placed on the outer surface of the test cylinder using clamps, the test cylinder 
being already prepared as it normally would for an experimental run. The test cylinder was 
then suspended horizontally and power was supplied to the cartridge heater. Once all four 
of the thermocouples reached a steady-state temperature, the temperatures were recorded 
along with the voltage and current being supplied to the heater. The equivalent thermal 
resistance was then derived using the equation 


R 


eq 


T - T 

C _ S 

IV 


(A.1) 


Temperatures which were to be comparable with those used during the experiments were 
obtained in a still air environment. Some additional tests were performed using a fan 
blowing air across the test cylinder to create a higher heat transfer rate away from the 
cylinder so that higher power levels to the heater could be reached while staying within the 
limits of the expected temperature range. An average temperature was then obtained using 
all three thermocouples, and this value was in turn utilized to derive an equivalent thermal 
resistance for each test run. The final value of the thermal resistance resulted in averaging 
the values from each test mn, obtaining 1.022 KAV with a maximum deviation of 0.089 
KAV. The results are annotated in Table 3. 

The same test results were used to determine the change in temperature along the 
length of the test cylinder. The difference in temperature along the length at the lower power 
settings demonstrated a fairly even distribution of heat, with a variation of less than 0.6°C. 
This variation increased as the power to the heater was also increased, as could be expected. 
For the temperature range in which the experiments were run (again, 26° to 40° C), a 
maximum variation of 1.6°C occurred, although this was at a very high level of power being 
supplied to the heater and was not representative of the remaining data. 
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Current 

(Amps) 

Voltage 

(Volts) 

Power 

(Watts) 

Interior 

Temp 

(C) 

Thermo¬ 
couple “A” 

(C) 

Thermo¬ 
couple “B” 

(C) 

Thermo¬ 
couple “C” 

(C) 

Average 

Surf Temp 

(C) 

Equiv 

Resistance 

(C) 

0.08 

8.30 

0.664 

38.4 

37.8 

38.1 

37.3 

37.73 

1.009 

0.07 

7.40 

0.518 

36.0 

35.5 

35.5 

35.4 

35.47 

1.023 

0.08 

8.40 

0.672 

37.3 

36.9 

36.8 

36.3 

36.67 

0.938 

0.09 

9.50 

0.855 

40.8 

40.2 

40.1 

39.5 

39.93 

1.018 

0.09 

10.30 

0.927 

42.9 

42.3 

42.0 

41.5 

41.93 

1.046 

0.10 

10.70 

1.070 

44.8 

44.2 

43.7 

43.2 

43.70 

1.028 

0.08 

9.00 

0.720 

27.4 

26.9 

26.6 

26.3 

26.60 

1.111 

0.09 

10.00 

0.900 

28.6 

27.9 

27.8 

27.2 

27.63 

1.078 

0.10 

11.00 

1.100 

30.1 

29.3 

29.2 

28.6 

29.03 

0.973 

0.11 

12.00 

1.320 

31.5 

30.5 

30.5 

29.6 

30.20 

0.985 

0.14 

15.50 

2.170 

37.6 

35.7 

36.0 

34.4 

35.37 

1.028 


Table 3: Thermal Resistance and Linear Temperature Distribution Trials 
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APPENDIX B. UNCERTAINTY ERROR ANALYSIS 


In order to obtain a measure of the reliability of the data obtained during 
experimentation, an uncertainty error analysis was performed for both the Nusselt number 
and the streaming Reynolds number. The analysis consisted of finding the maximum 
possible deviation for all components appearing in the equations which define both 
parameters. Then, in a standard fashion, a root mean square analysis was performed to 
derive a reasonable overall possible error for each experimental run. The error analysis 
formulae were themselves incorporated into the “results worksheet” provided in Appendix 
C so that each data point has an associated possible error derived from the input provided. 

The Nusselt number can be derived as shown from Eqs. 37 through 40 


Nu 


_TY_ 

Tikl(T - IVR - T) 

^ c eq a'' 


(B.l) 


where I is the length of the test cylinder. 

Using the manufacturers’ recommended equipment error ranges, the maximum 
uncertainty in each measured component in Eq. B.l is as follows 


V = 0.05% reading + 0.02% full scale + 1 digit 

(B.2) 

I = 0.5% RDG + 1 digit 

(B.3) 

r, = ±0.5" c 

(B.4) 

= ±0.5" C 

(B.5) 

The maximum uncertainty in the equivalent resistance is 

calibration data in Appendix A as 

obtained from the 
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R = ±0.089 K/W 

eq 


(B.6) 


The analysis for the overall uncertainty itself is structured as follows. For the Nusselt 
number, the root mean square error is given by 


a(A^m) = 

‘ aA; 


(B.7) 


where X; represents each individual component of Eq. B.l. The-term, as the partial 

derivative indicates, physically represents the sensitivity of the Nusselt number to tiie 
variable Xj, provided all other variables are unchanged. The aX. represents the uncertainty 
in the corresponding variable as given in Eqs. B.2 to B.5. For instance, the contribution to 
the uncertainty due to the voltage measurement is 




^^aV = Nu —— 

ay y(r. - - r) 


(B-8) 


A value with more significance, though, is the individual fractional uncertainty which 
takes into account the calculated value of the Nusselt number and can be expressed as a 
percentage possible error. This is represented by dividing the individual uncertainty by the 
Nusselt number in the following manner 


individual fractional uncertainty = - = - 

Nu V(T^ - VIR^^ - TJ 


(B.9) 


This new term leads directly to the desired method of expressing the possible error 
in the calculated value of the Nusselt number as an overall fractional uncertainty using a root 
mean square analysis. 


a(Nu) _ V ^ 

Nu i Nu 


(B.IO) 


In a similar fashion, the individual fractional uncertainty of the remaining terms are 
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as follows 


(aNu), _ (T^. - TJiAP) 


Nu 

I(T, 

- yiK, - 

TJ 

(aNu).j. 




Nu 


- ■ 


(aNu)j. 


(aT) 


Nu 

(r 

- y^K - 


(^aNu)^ 




Nu 

(Tc 

- yiRe, - 

- TJ 


(B.ll) 

(B.12) 

(B.13) 

(B.14) 


By examining the data results, the largest contributor to the overall error in the Nusselt 
number is due to the current. This is caused by the small currents being utilized during the 
experiment, which were as low as 0.05 amps. Since the uncertainty in the current is of the 
order of 0.01 amps, there can be an error in the calculated Nusselt number of approximately 
20% due to the current term alone. One way to lessen the effect that the current term has on 
the overall error is by decreasing the voltage output of the power supply and hence increasing 
the current needed to maintain the same power being generated. This capability, though, is 
not a feature of the equipment being used. It must also be noted that although 
= ±0.5" Cin Eqs. B.4 and B.5, the calibration of the thermocouples described in 
Appendix A indicated an error of less than 0. TC for the temperature range of the experiment 
If this is taken into account, the error due to the ambient temperature (TJ and the center 
temperature (Tj.) terms in Eq. B.l would diminish by 80%. 

The error analysis for the streaming Reynolds number is similar to that just performed 
for the Nusselt number. Utilizing Eq. 17 and substituting Eqs. 33 and 41 into it, the value 
of Rj can be shown to be 


RjT^ + 

27l^YP„'(0.523)2 


(B.15) 


Recall that VJsJl is the multimeter output voltage as derived from the pressure transducer 


after being passed through the 100 gain preamplifier. Thus, Rj becomes a function of only 
the following measured variables: the frequency, the ambient temperature and the pressure 
transducer output as read on the multimeter, while all additional parameters remain constant. 
The maximum uncertainty for the thermocouple is the same as previously listed (±0.5°C) 
while the other two variables have uncertainties provided in the manufacturers’ 
specifications. Calibration data for the pressure transducer used during the experiment 
indicates an uncertainty equivalent to 0.15% of the percentage of Full Scale Output (FSO) 
where the FSO is 254 mV. This correlates to a maximum deviation of 2.3 mV after passing 
through the 100 gain preamplifier. In addition to this, the multimeter which is used to read 
this voltage has an error of 0.06% reading + 0.03% Range. These combine to give a total 
microphone voltage output error of 

Mic = 2.3 mV + 0.06% reading + 0.03% range (B.16) 

The uncertainty in the frequency signal from the function generator error is given as 20 ppm, 
i.e., 

/ = 20x10^^ reading (B.17) 


Again a root mean square analysis was performed in a manner similar to Eq. B.7 to derive 
the overall fractional uncertainty for the streaming Reynolds number. 


^ - rE 
R i R. 


(B.18) 


The individual fractional uncertainties are 

R^ (T^ + 273.15) 
^ 2{AMic) 


(B.19) 

(B.20) 
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(B.21) 


Rs f 

The largest contribution to the error in the streaming Reynolds number is the error 
due to the microphone output voltage, specifically, the possible error in the pressure 
transducer itself. However, the data presented in Appendix C shows that this error is very 
small and in the range of <2% of the total value. 
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APPENDIX C. EXPERIMENTAL DATA 


The following pages contain the data obtained through experimentation in 
spreadsheet fashion. All of the relevant parameters are listed, although some constants have 
been left out due to size constraints. 
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